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Abstract:  When  data  comes  from  a  distribution  belonging  to  the  domain  of 
attraction  of  a  stable  law,  Athreya  (1987a)  showed  that  the  bootstrapped  sample 
mean  has  a  random  limiting  distribution,  implying  that  the  naive  bootstrap 
could  fail  in  the  heavy-tailed  case.  The  goal  here  is  to  classify  all  possible 
limiting  distributions  of  the  bootstrapped  sample  mean  when  the  sample  comes 
'rom  a  distribution  with  infinite  variance,  allowing  the  broadest  possible 
wetting  for  the  (nonrandom)  scaling,  the  resample  size,  and  the  mode  of 
convergence  (in  law).  The  limiting  distributions  turn  out  tc  '~>e  infinitely 
divisible  with  possibly  random  Levy  measure,  depending  on  the  3sample  size. 

An  averaged-bootstrap  algorithm  is  then  introduced  which  eliminates  any 
randomness  in  the  limiting  distribution.  Finally,  it  is  shown  that  (on  the 
average)  the  limiting  distribution  of  the  bootstrapped  sample  mean  is  stable  if 
and  only  if  the  sample  is  taken  from  a  distribution  in  the  domain  of  (partial) 
attraction  of  a  stable  law. 
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1.  INTRODUCTION 

The  bootstrap  (Efron  (1979))  is  a  resampling  technique  for  nonpara- 

metrically  modeling  the  distribution  of  a  general  statistic  T  =  T  (ft  ; F) , 

n  nv  n 

where  ft^  =  (Xj,...,X  )  is  a  random  sample  of  size  n  from  a  population 

distribution  F.  Bootstrap  "replicates"  T*  =  T  (3/  ;F  )  of  the  statistic  are 

n  n 

generated  by  resampling  *?/  =  (Y,,...,Y  )  from  the  empirical  distribution  F  of 

m  1  m  '  n 

n  n 

(Xj,...,Xn),  i.e.,  the  Y^’s  are  conditionally  iid  with  distribution  having  mass 

—  at  each  of  the  points  X, , . . .  ,X  .  The  conditional  distribution  S£(T  Ift  ), 
n  1  n  v  n1  ny 

which  is  necessarily  random ,  is  intended  to  approximate  the  true  nonrandom 
distribution  S£(Tn). 

In  the  present  paper  we  deal  only  with  the  simplest  statistic  T  ,  the 
standardized  sample  mean.  The  sample  mean  is  the  easiest  statistic  to  handle 
analytically,  and  if  the  bootstrap  algorithm  fails  for  this  simple  statistic, 
we  cannot  expect  it  to  succeed  for  more  complex  statistics.  Our  work  is 
motivated  by  the  results  of  Bickel  and  Freedman  (1981),  Singh  (1981),  Athreya 
(1986a, b,  1987a, b),  Gine  and  Zinn  (1990a, b),  Arcones  and  Gine  (1989),  Hall 
(1990),  and  Knight  (1989). 

We  denote  by  H  (x)  the  conditional  distribution  function  of  T*  given  ft  , 
nv  *  n  &  n 

H  (x)  =  P(T  £  x  |  ft  ),  where  T’  =  a  2,  ,(Y,-c  )  is  the  bootstrap  version  of 
n  n  n'  n  n  i=lv  in 

the  sample  mean  with  possibly  data-driven  scaling  an  and  centering  c^  (cn  =  X^ 
when  F  has  finite  mean,  and  otherwise  cn  =  0),  and  with  (nonrandom)  resample 
size  mn  (which  may  differ  from  the  original  sample  size  n).  As  Hn(x)  is  a 
random  distribution  function,  several  modes  of  convergence  are  possible  (such 
as  convergence  a.s.,  in  probability,  in  law),  and  the  limit  may  be  random  or 
nonrandom .  We  judge  the  performance  of  the  bootstrap  by  asymptotically 
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comparing  H  to  S£ ,  the  limiting  distribution  of  the  standardized  sample  mean 
T  .  We  now  review  some  of  the  known  results  on  the  bootstrap.  Each  of  these 
results  involves  specific  assumptions  on  an  and  cn>  for  which  the  reader  is 
referred  to  the  original  papers. 

A.  (Sufficient  conditions.! 

(A.l)  (Bickel  and  Freedman  (1981))  If  F  has  finite  variance  and  m^  — * 09 . 

then 


sup  |H  (x)  -  5£w(x)  |  g£*  0; 

x€R  n 

here  is  the  standard  normal  distribution. 

00 

(A. 2)  (Arcones  and  Gin6  (1989))  If  F  belongs  to  the  domain  of  attraction 

of  the  normal  law  and  m  — *  «,  then 

n 

sup  |H  (x)  -  (x)  |  0; 

X€P.  n  00  nH“ 

here  S£w  is  the  standard  normal  distribution.  (See  Athreya  (1987b)  for  mn=n.) 

(A. 3)  If  F  belongs  to  the  domain  of  attraction  of  an  a-s table  law,  0<a<2, 

then 

(i)  (Athreya  (1987a))  if  n>n=n,  the  bootstrap  fails  (in  fact  Hn(x) 
converges  in  law  to  a  random  distribution  which  differs  from  $m,  an  a-stable 
distribution) ; 

(ii)  (Athreya  (1987b))  if  mn/n  — »  0,  the  bootstrap  works  in  probability, 
i.e.,  the  conclusion  of  (A. 2)  holds  with  an  a-stable  law  $ L; 

(iii)  (Arcones  and  Gine  (1989))  if  (mn/n)  loglog  n  — >  0.  the  bootstrap 
works  a.s.,  i.e.,  the  conclusion  of  (A.l)  holds  with  an  a-stable  law  J£w. 

(A. 4)  (Hall  (1990))  If  F  has  slowly  varying  tails  and  one  tail  completely 
dominates  the  other,  and  mn=n,  then  the  bootstrap  fails  (in  fact  Hn  converges 
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in  probability  to  a  nonrandom  Poisson  distribution). 

B.  (Necessary  conditions . ) 

(B.l)  (Arcones  and  Gin6  (1989))  If  inf  mn/n  >  0  and  for  some  possibly 
random  distribution  H(x) , 

H  (•)  H(»)  a.s. 

then  F  has  finite  variance,  and  H  is  S£w  (the  standard  normal  distribution). 

(See  Gine  ana  Zinn  (1990b)  for  n>n=n) . 

(B.2)  (Arcones  and  Gin6  (1989))  If  mR  — *  an  are  nonrandom,  and  for 

some  possibly  random  distribution  H(x) , 

Hn(0^H(-)  in  probability. 

then  F  belongs  to  some  domain  of  partial  attraction  and  H  is  (an  infinitely 

divisible  distribution).  If  furthermore  lim  inf  ^  m  /n  >  0,  then  F  belongs  to 

n-#°  n 

the  domain  of  attraction  of  a  normal  law,  and  H  is  (the  standard  normal 
distribution) . 

(B.3)  (Hall  (1990))  If  n>n=n,  aR  are  possibly  random,  and  for  some 
nonrandom  distribution  H(x), 

H  (•)  — *.H(*)  in  probability  (at  continuity  points  of  H) , 
n  n-Hoo 

then  either  F  belongs  to  the  domain  of  attraction  of  a  normal  law  and  H  is 
(the  standard  normal  distribution),  or  F  has  slowly  varying  tails  with  one  tail 
completely  dominating  the  other  and  H  is  a  Poisson  distribution. 

Our  interest  is  to  extract  information  from  in  the  broadest  possible 
setting,  i.e.,  with  minimal  restrictions  on  the  underlying  distribution  F,  the 
resample  size  mn>  and  the  standardization  a^  and  c^,  and  allowing  the  weakest 
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form  of  convergence  (in  law)  to  a  possibly  random  limit.  Our  goal  is  to 
classify  all  possible  limiting  distributions,  and  then  to  eliminate  any 
randomness  in  the  limit. 

When  F  is  general  with  infinite  variance,- assuming  nonrandom  an>  we  give 
in  Section  2  necessary  and  sufficient  conditions  for  weak  convergence  of  the 
finite  dimensional  distributions  of  (Hn(x),  -®<x<«),  and  we  classify  all 
possible  limits.  The  result  (Theorem  1)  shows  that  Hn  may  have  a  limit  even 
when  the  sample  mean  does  not,  and  that  the  limit  of  Hn  is  always  infinitely 
divisible  with  possibly  random  Levy  measure,  depending  on  the  choice  of  the 
resample  size  mn> 

For  the  case  where  H  (x)  has  a  random  limit,  we  introduce  in  Section  3  an 
averaged-bootstrap  algorithm  and  show  in  Theorem  2  that  the  averaged-bootstrap 
distribution  converges  a.s.  to  a  nonrandom  limit,  G(x),  which  is  the  expected 
value  of  the  random  limit  of  Hn(x).  If  G  contains  useful  information  about 
then  averaged-bootstrap  resampling  is  useful,  i.e.,  bootstrap  resampling  is 
still  worthwhile  on  the  average. 

In  Section  4  we  apply  these  results  to  population  distributions  in  the 
domain  of  attraction  of  a  stable  law.  We  show  that  the  random  limiting 
distribution  of  Hn(x)  obtained  by  Athreya  (1987a)  with  a  data-based  scaling,  is 
not  stable  even  on  the  average  (Theorem  4).  However,  with  a  nonrandom  scaling, 
the  averaged  bootstrap  actually  yields  a  stable  distribution  with  the  same 
index  as  but  with  possibly  different  skewness  and  scale  parameters  (Theorem 
5).  Appropriately  adjusted,  the  averaged-bootstrap  algorithm  works  a.s.  in 
some  cases,  i.e.,  yields  the  correct  skewness  and  scale  of  (Theorems  6  and 
7).  Moreover  the  limit  distribution  G  of  the  averaged  bootstrap  is  stable  if 
and  only  if  F  belongs  to  the  domain  of  (partial)  attraction  of  a  stable  law 
(Theorem  8). 
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2.  GENERAL  POPULATION  DISTRIBUTION  WITH  INFINITE  VARIANCE 

A  fundamental  question  is:  What  are  the  possible  limits  of  the 
bootstrapped  conditional  distribution  ^(T*|£n)  when  we  do  not  have  any 
restriction  on  the  population  distribution  F,  the  resample  size  mn»  and  the 
scaling  aR.  We  consider  on  the  one  hand  general  population  distributions  F 
with  infinite  variance  (i.e.,  without  the  domain  of  attraction  restriction  of 
results  A),  and  on  the  other  hand  we  allow  general  resample  size  mR  and 
convergence  in  distribution  of  Hn  (rather  than  the  more  restrictive  a.s.  or  in 
probability  convergence  of  results  B). 

The  following  theorem  gives  necessary  and  sufficient  conditions, 
explicitly  in  terms  of  the  univariate  distribution  F  of  X,  for  the  convergence 
in  law  of  S£(T*|iXn)  and  describes  the  possible  limits. 

Theorem  1 .  Assume  the  population  distribution  F  has  infinite  variance,  a^  is 
nonrandom , 

_o 

(1.0)  m  — ►  »,  a  — >  00 ,  m  a  — »  0  as  n  — ►  «, 

v  '  n  n  n  n 

and  when  F  has  finite  mean,  assume  furthermore 


(1.1)  lim  lim  n  P{|x|  >  a  max(M,en/m  )}  =  0  for  all  e  >  0, 

Mh“  n-*o  n  n 


(1.2) 

Then 

(1.3) 


lim  lim  m  a  |x|l( |xj  >  Ma  ,  |x|  <  ea  n/m  )}  =  0  for  some  e  >  0. 
M-*o  n-*°  n  n  n  n  n 


{Hn(Xl) . H^)}  (H(x1),...,H(xk)} 


for  some  possibly  random  distribution  function  H(x),  any  integer  k  and 

Xj . x^  €  IR,  if  and  only  if  the  following  three  conditions  (l.a)-(l.c)  are 

satisfied: 

(l.a)  There  exists  a  Levy  measure  X,  such  that  for  all  y>0, 
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for  every  6  >  0  such  that  v{6)  =  0. 

.(l.c)  For  all  6  >  0  such  that  u{6)  =  0,  there  exists  a ^  such  that 

lim  lim  mnan2  E{X2I[|X|  £  anmin( e , Vfin/m^) ] }  =  a2. 
e-O  n-i» 

In  this  case,  the  possibly  random  characteristic  function  <f^( • )  corresponding 
to  H(*)  is 

exp{-  |t 2ct2  +  /^[e1  tx-l-i tT(x) ]dX(x)  }  if  — *  0, 

(1.4)  <PH(t)  =  ‘  exp{-  |-t2a2  +  «£,  [eitx-l-itT(x)]dN(x)  }  if  —•  — »  c  €  (0,«), 
exp{-  ^t2a2  -  jjj^W  }  otherwise, 

2  2 

where  a  =  lim^Q  a g,  t(x)  =  x  if  F  has  finite  mean,  and  r(x)  =  0  otherwise,  N 

-1  2 

is  a  Poisson  random  measure  with  intensity  measure  c  X,  a  +W  £  0,  and  the 
characteristic  function  <p^(u)  of  W  is  expl/Q^^-lJduCx)} .  Furthermore,  the 
expected  value  G(x)  of  H(x)  has  characteristic  function  <#>g(t)  (=  Ef>^(t))  which 
is  always  infinitely  divisible  with  the  following  form 
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(1.5)  exp{-  \t2a2  +  J'^)[eitX-l-itT(x)]dX(x)  } 


exp{-  it 22  +  c-1/” {exp[c(eitX-l-itT(x))] 


m 

lfr-0' 

m 

l}dX(x)}  if  c€(0.«) 


exp{-  it2<72}  Eexp{-  it2!} 


otherwise. 


Theorem  1  classifies  all  possible  weak  limits  of  ^(T*l®n)-  They  are 

either  the  law  of  an  independent  sum  of  a  Gaussian  r.v.  and  a  compound  Poisson 

r.v.  with  possibly  random  Levy  measure,  or  else  the  law  of  an  infinitely 

divisible  scale  mixture  of  a  Gaussian  r.v.,  i.e.,  a  Gaussian  law  with  zero  mean 

and  random  infinitely  divisible  variance.  The  (possibly  random)  variance  of 

the  Gaussian  component  and  the  L6vy  measure  of  the  compound  Poisson  component 

depend  on  the  choice  of  resample  size  and  standardization. 

Conditions  (1.1)  -  (1.2)  ensure  that  the  tails  of  the  integrals  in  (1.4) 

when  F  has  finite  mean  are  eventually  negligible.  The  Gaussian  component  of 

2 

the  limit  has  variance  a  determined  by  the  asymptotic  behavior  of  a  certain 
truncated  second  moment  of  F,  described  in  (l.c).  The  compound  Poisson 
component  of  the  limit  is  determined  by  the  tail  behavior  of  F  described  in 
(l.a).  And  the  Gaussian  mixture  has  an  infinitely  divisible  variance  whose 
Levy  measure  v  is  determined  by  the  asymptotic  behavior  of  a  certain  truncated 
probability  described  in  (l.b). 

The  form  of  the  limits  (1.4)  and  (1.5)  can  be  interpreted  as  follows.  The 

bootstrap  resample  comes  from  the  empirical  distribution  F  .  If  we  take  the 

resample  size  mn  much  smaller  than  the  original  sample  size  n  (mn/n  — *  0)»  then 

from  the  perspective  of  the  resample  V  (mn  large),  F^  is  essentially  like  the 

n 

nonrandom  population  distribution  F.  Thus  the  distribution  of  the  bootstrapped 
sample  mean  will  yield  a  deterministic  limit  if  it  exists  (as  in  (1.4)).  This 
phenomenon  appears  in  the  literature  in  other  contexts  (see  Bretagnolle  (1983) 
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and  Swanepoel  (1986)). 

Another  Interesting  effect  arises  from  the  fact  that  Fr  has  finite  support 

(trivially  symmetric  tails).  If  we  take  the  resample  size  mR  much  larger  than 

the  original  sample  size  n  (m  /n  — *  m) ,  then  V  is  like  a  very  large  sample 

n  m 

n 

from  Fn  (rather  than  F).  Hence  the  structure  of  Fr  (symmetric  tails)  will  be 
reflected  in  the  limiting  distribution  H,  which  indeed  is  always  symmetric  in 
this  case  (provided  it  exists). 

Under  conditions  (l.a)-(l.c),  the  limit  H,  or  its  expected  value  G,  may 

contain  some  information  about  This  suggests  that  Hn>  or  its  expected 

value  EH  ,  may  contain  some  information  about  i£(T  ).  In  Section  3.  we  develop 
n  n 

the  averaged-bootstrap  algorithm  in  order  to  extract  this  information. 

There  are  many  distribution  functions  with  infinite  variance  whiA  do  not 

belong  to  the  domain  of  attraction  of  any  stable  law,  i.e.,  whose  sample  mean 

does  not  have  a  limiting  distribution.  Theorem  1  suggests  the  existence  of  a 

limit  for  S£(T*|9t  )  in  some  such  cases. 
v  n‘  n' 

3.  AVERAGED  BOOTSTRAP 

In  Section  2  we  have  classified  all  possible  limiting  (cond  ional) 

distributions  of  T  .  It  is  clear  from  Theorem  1  that  the  limits  may  be  random, 

and  therefore  different  from  the  nonrandom  limit  if  of  S£(T  )  when  it  exists, 

co  v  n7 

e.g.,  when  F  belongs  to  some  stable  domain  of  attraction.  The  question  arises: 

How  can  we  extract  useful  information  about  %m?  For  this  purpose  we  now 

introduce  the  averaged-bootstrap  algorithm. 

Partition  the  data  set  into  2  blocks,  each  having  k  observations, 

n  n 


,Xk  ,Xk  +1’ 
n  n 


,X2k . X(2  -l)k  +1’ 

n  v  n  J  n 


,X«k 
n  n 


*k 


n 


2^k 


n 


2  \ 
n  n 


n=£  k  ,  as  follows: 
n  n 
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*  -1  n 

Instead  of  bootstrapping  the  entire  sample  to  obtain  Tn=  aR  2^  (Y^-c^) ,  we 
apply  the  bootstrap  algorithm  within  each  block  to  obtain  a 

block-bootstrapped  version  of  the  sample  mean 


-1  „  n 


=  J\  Ji=l  ‘jYi  '  A,’’ 


th 


where  {jYj>  }  are  resampled  from  the  j  block  ,  i.e.,  they 

conditionally  independent  with  common  distribution 

Vi  * xu-i)k+i  I  j\>  ■  r  •  «*»>,• 


are 


is 


The  conditional  distribution  of  the  block-bootstrapped  sample  mean  .T,  i 

J  n 

denoted  by 

Anw  ■  Vkn  i  *  l  jV 

In  T  we  use  c  =  X  if  F  has  finite  mean,  c  =  0  otherwise;  and  in  .T.  we  use 
n  n  n  n  j  k 


=  0 


j°k  =j^k  sample  mean  in  the  j  block)  if  F  has  finite  mean,  ^c^ 

otherwise.  For  a  data-driven  scaling,  we  use  ^a^  =  a^  (^3^  )  in  , 

analogously  to  using  a  =  a  (9t  )  in  T  . 

n  nv  n'  n 

Note  that  {  .H,  (x),  l<j£5  }  are  iid  random  distributions;  their  average  is 
j  Kn  n 


denoted  by 


-  f  j"k  <*>’ 

n  j=l  n 


and  is  referred  to  as  the  averaged-bootstrap  distribution. 

We  now  show  that  the  averaged-bootstrap  algorithm  converges  a.s.  when  the 
standard  bootstrap  algorithm  converges  in  law  (the  weakest  possible  mode). 
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Theorem  2.  If  l  =  n  .  0<6<1.  and  for  each  x  €  IR, 
-  n 


H(x)  •—+  H(x) , 

n-*» 


where  H(x)  is  a  possibly  random  distribution  function,  then  for  each  x  €  !R 


H  (x)  ^JL_*G(x), 
nv  J  n-*»  v  J 


where  G(x)  =  EH(x).  Furthermore  if  G(x)  is  continuous,  then  the  convergence  is 
uniform  in  x. 


Thus  the  averaged-bootstrap  algorithm  produces  a  nonrandom  limiting 

distribution  G(x),  and  it  is  of  interest  to  investigate  the  relationship 

between  G  and  i£m.  If  G  and  i£m  agree,  then  the  averaged  bootstrap  works  a.s. 

and  H  contains  useful  information  about  S£(T  ). 
n  v  nJ 

4.  APPLICATIONS  TO  POPULATION  DISTRIBUTIONS  IN  DOMAINS  OF  ATTRACTION 

In  this  section  we  apply  Theorems  1  and  2  to  population  distributions  in 
the  domain  of  attraction  of  a  stable  law.  Each  of  the  stated  results  (except 
Theorem  3)  includes  two  parts.  First  we  establish  convergence  of  the  random 
distribution  Hn(x)  to  a  possibly  random  limiting  distribution  H(x),  using 
Theorem  1  when  appropriate.  Secondly,  an  explicit  expression  for  G(x)  =  EH(x) 
is  obtained,  again  using  Theorem  1.  By  comparing  G  with  i£m,  we  can  in  each 
case  determine  whether  the  averaged  bootstrap  (Hn)  is  of  any  practical  use  (via 
Theorem  2). 

Before  stating  our  results  we  briefly  review  the  stable  laws  S(a,/J,s)  with 
index  a  €  (0,2],  skewness  P  €  [-1.1].  scale  parameter  s  £  0  (the  actual  scale 

1/tt 

is  sA/  ;  here  for  simplicity  we  use  s  instead),  and  characteristic  function: 
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*  exp{-s|t|a[l  -  iptan(a7r/2)sgn(t)]  },  a  j*  1, 
exp{-s | t |[1  +  ip(|)log|tisgn(t)]  },  a  =  1. 


A  distribution  function  F  belongs  to  the  domain  of  attraction  of  a  stable  law 

with  0<a<2,  denoted  F  €  D(S(a,/3,s)) ,  if  and  only  if  there  exist  constants 

Cj.Cg  £  0,  and  a  slowly  varying  function  L(x)  such  that  l-F(x)  ~  c^x  ^(x), 

F(-x)  ~  CgX  ^(x),  as  x  — ►  »,  where  (3  =  (c^-CgJ/fcj+Cg)  •  F  belongs  to  the 

domain  of  attraction  of  a  Gaussian  law  (a=2)  if  and  only  if  E[X^I( jx|^x)]  is  a 

slowly  varying  function  L(x).  F  belongs  to  the  normal  domain  of  attraction  of 

a  stable  law,  denoted  by  F  €  D^(S(a,/3,s)) ,  if  and  only  if  F  €  D(S(a,/3,s))  with 

L(x)  =  c  for  some  constant  c  >  0.  These  are  precisely  the  distribution 

functions  whose  sample  means  (properly  normalized)  have  limiting  distributions, 

i.e.,  F  €  D(S(a,p,s))  (respectively  D^(S(a,p,s)))  if  and  only  if  there  exist 

>  0  and  C  €  IR  such  that  A  *(nX  -C  )  converges  in  law,  in  which  case  the 

limiting  distribution  S£  is  S(a,/3,s)  and  A  =  n^^fn)  (respectively 

00  n 

*  i/a\ 

A  =  n  ). 

When  F  €  D(S(a,p,s)),  it  can  be  easily  verified  that  the  statements  (A. 2) 
and  (A. 3)  (ii)  in  Theorem  A  can  be  derived  from  Theorem  1.  Namely  there  exist 
appropriate  scaling  a^  and  resample  size  mn  such  that  the  bootstrapped  sample 
mean  has  asymptotic  limiting  distribution  as  given  in  Theorem  1,  depending  on 
the  resample  size.  On  the  other  hand,  if  the  bootstrapped  sample  mean  has  a 
limiting  distribution,  then  the  limit,  the  scaling  and  the  resample  size 
necessarily  must  satisfy  certain  conditions  described  as  follows. 


Theorem  3.  Assume  F  €  D(S(a,p,s))  has  infinite  variance  and  assumptions 
(1.0)-(1.2),  (1.3)  hold.  Then  (an)  satisfies 

lim  m  a  a  L(a  )  =  constant  if  m  /n  — »  c  €  [0,“), 

_  n  n  v  n''  n 
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lim  n(a  v4i/m  )  “Lfa  vii/m  )  =  constant  if  m  /n  — ► 00 . 

^  v  n  ny  v  n  n'  n 

n-J® 

Moreover,  if  a=2,  then  H  is  a  nonrandom  normal  distribution. 

Now  we  examine  the  performance  of  the  averaged  bootstrap  with  Athreya’s 
(1987a)  data-driven  scaling. 

Theorem  4.  If  F  €  D(S(a,p,s)),  a  /  2,  mR  =  n,  and  Athreya’s  data-driven 

scaling  a  (St  )  =  X  (the  sample  maximum)  is  used,  then  (1.3)  holds  and  G(x) 
n  n  n  »  n 

has  characteristic  function 

Cj  exp(f t(l)) 

V>G(t)  =  - ; -  . 

C1  +  Cat1  "  exp(ft(x))]dAa(x) 

where  ft(x)  =  e-tx-l-ita(x) ,  a(x)  =  r(x)  if  a  /  1,  a(x)  =  t(x)I(|x|<1)  if 

a  =  1,  and  the  measure  X  satisfies  for  all  v.  N  0, 

a 

X  [x,00)  =  c.x  a,  X  '(-«,-x]  -  cnx  a. 
ar  '  1  av  J  2 

Since  is  quite  different  from  the  characteristic  function  of  i£m,  the 
bootstrap  fails  even  on  the  average  with  this  scaling.  Note  that  this 
data-based  scaling  results  in  a  bootstrap  algorithm  which  uses  no  specific 
features  of  the  population  distribution  F  other  than  the  knowledge  that  it 
belongs  to  the  domain  of  attraction  of  some  stable  law.  However  with  an 
appropriate  nonrandom  scaling,  the  averaged-bootstrap  algorithm  can  be  improved 
in  the  sense  that  G  is  stable  with  the  same  index  as  i£m  but  with  generally 
different  skewness  and  scale. 

Theorem  5.  If  F  €  D(S(a,/3,s)) ,  1  /  a  €  (0,2),  mn  =  n  and  an  is  a  nonrandom 
scaling  based  on  the  population  distribution  as  follows: 

n  P(X  £  a  )  — ♦ 

y  -  nJ  n-*» 


1. 
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then  (1.3)  holds  and  G(x)  is  Sfa.pdgCct) ,sd.^(a)) ,  where 


r  e-1 


dx(“)  =  \ 


«  ,  a 

y  k  +  i 

Lk=i  <k*ni  ] 

,  l<a<2. 

«  ,  a 

2  — 
k'  ’ 
k=l 

0<a<l, 

1  _  2 
ed j(a)  ’ 

l<a<2. 

1. 

0<a<l . 

When  o=l,  G(x)  is  not  Cauchy  even  with  the  above  nonrandom  scaling,  and 
its  characteristic  function  has  a  form  more  complex  than  that  in  Theorem  3  and 
is  therefore  not  shown  here. 

Note  that  in  each  of  the  intervals  (0,1)  and  (1,2),  d^(a)  is  increasing, 
infinitely  differentiable  with  all  derivatives  positive  and  d^(0+)  =  1-e 
dj(l-)  =  1,  dj(l+)  =  2/e,  d^(2-)  •■=  1;  see  Figure  1  (c=l).  Therefore  the  scale 
of  G  is  smaller  than  that  of  The  factor  dgfa)  is  also  increasing  on  (1,2) 

with  d2(l+)=0,  =  1-2/e;  see  Figure  2  (c=l).  Hence  the  skewness  pdgfo) 

of  G  is  the  same  as  the  skewness  p  of  t£m  for  a  €  (0,1),  and  has  the  same  sign 
but  is  reduced  by  a  factor  of  at  least  .736  for  a  €  (1,2). 

When  F  belongs  to  the  normal  domain  of  attraction  of  a  stable  law,  the 
bootstrap  algorithm  can  in  some  cases  be  appropriately  modified  so  that  the 
averaged  bootstrap  will  work.  The  simplest  cases  are  when  0<a<l  (so  there  is 
no  reduction  in  skewness)  or  when  the  skewness  is  known,  e.g.,  in  the  symmetric 
case.  Then  only  a  scale  adjustment  to  the  algorithm,  depending  on  the  index  a, 
is  necessary;  when  a  is  unknown,  an  appropriate  estimate  may  be  used. 

Theorem  6.  If  F  €  Djj(S(a,/3,  s)) ,  mn=n,  either  0<a<lor/3  =  0  and 

i/a  (a  ) 

1  /  a  €  (0,2),  sind  an(®n)  =  [n  dj(an(SCn))]  ,  where  an(&n)  is  311  estimate 


M 


A  P 

of  a  satisfying  (an~a)  log  n  >  0,  then  (1.3)  holds  and  G  =  both  being 
S(a,p,s) . 


For  instance  the  estimate  of  a  given  in  Hall  (1982)  works  in  Theorem  6. 

Notice  the  skewness  of  G  in  Theorem  5,  Pq  =  pd2(ct),  is  always  smaller  in 
absolute  value  than  the  skewness  P  of  when  l<a<2.  Hence  if  prior  to 
bootstrapping  we  transform  the  original  data  into  new  data  belonging  to  the 
domain  of  attraction  of  a  stable  law  with  skewness  pd2*(a),  and  if  we  then 
appropriately  adjust  the  standardization  as  in  Theorem  6,  using  estimates  of  a 
and  p  if  necessary,  we  obtain  a  bootstrap  algorithm  which  works  on.  the  average 
provided  |p 1  £  «.  i) . 

A  A 

Theorem  7.  Assume  F  €  D^(S(a,p,s)) ,  1  ^  a  €  (0,2),  |p|  £  d2(a),  and  an>Pn  are 
consistent  estimates  of  a  and  p,  satisfying 

(a  -a)  log  n  — 0  if  0  <  a  <  1, 
v  n  J  n-*° 

l-l/a,*  \  P  n  1  -l/a,£  P .  r,  .j.  i  ✓  s  n 

n  (a  -a)  — 0,  n  (P  -p)  — -*  0  if  1  <  a  <  2. 
v  n  J  n-*»  Vfn  n-*» 


Transform  the  data  9C  to  Z  as  follows: 

n  n 


A  A 


Z.  =  b (a  ,P  )  XT  -  b  (a  ,p  )  X. 
i  v  n  rn'  i  '  n  'll'  i 


where 


b+  =  b+(a,p)  =  (l-p)[l+P/d2(a)],  b"  =  b  (a.p)  =  (l+P)[l-p/d2(a)]  , 

^  __  ^  _n 

and  consider  the  bootstrapped  sample  mean  T  =  a  2v  i(Y.  -  c  ),  where 

n  n  1=1  i  n 


an  ■  K(in)d-^)]1/an.  cn  =  * 


Z  if  a  >  1, 
n 

0  if  a  <  1, 


mn  =  n,  and  are  resampled  from  {Zj}^_j.  Then  (1.3)  holds  and  G  = 
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both  being  S(a,|3,s). 

For  instance  the  estimates  of  a  and  p  given  by  Zolotarev  (1986)  work  here. 

Up  to  now  we  have  discussed  in  detail  that  when  F  €  D(S(a,/3,s))  and 
n>n  =  n,  the  expected  value  G  of  the  limiting  distribution  H  is  a-stable,  with 
possibly  different  skewness  and  scale  from  and  how  to  make  bootstrap 
resampling  work  a.s.  by  using  the  averaged-bootstrap  algorithm.  The  following 
theorem  shows,  without  any  restriction  on  the  resample  size  mn,  that  G  is 
a-stable  if  and  only  if  F  belongs  to  the  domain  of  (partial)  attraction  of  an 
a-stable  law. 


Theorem  8.  Assume  (1.0)-(1.2)  and  (1.3)  hold. 

(i)  If  n>n/n  — ►  c  €  (0,09],  then  G(x)  is  a-stable  with  1  j*  a  €  (0,2]  if  and 
only  if  F  €  D(S(a,0,s)). 

(ii)  If  m  /n  — *  0,  then  G(x)  is  a-stable  with  1  ^  a  €  (0,2]  if  and  only  if 

F  belongs  to  the  domain  of  (mn)-partial  attraction  of  an  a-stable  law 

(F  €  (mn)-D  (S(a,p,s)) ,  i.e.,  A  *(m  X  -C  )  converges  in  law  for  some  A  >0 
n  P  n  n  mn  n  n 

and  €  R) . 

In  either  case  G(x)  is  Gaussian  with  mean  zero  and  variance  d(2,c)  if 
a  =  2,  otherwise 


G(x)  is 

where 


S(a,  p,  sd(a,0))  if  m  /n  — *  0, 

S(a,  /Jdgfa.c),  sd(a,c)dj(a,c))  if  mn/n — *c  €  (O*00)* 

S(a,  0,  sd(a,»)[2a^cos(a7r/4)]  *)  if  mn/n  — *  09 . 

lim  m  a  ^(a  )  if  m  /n  — >  c  €  [O,09), 

^nnvn'  n  Ly 

nH“ 

lim  n(a  'Jn/m  )  Lfa  vh/m  )  if  m  /n  — * 

_  v  n  ny  v  n  ny  n 

n-*° 


d(a,c)  = 


00 
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d^a.c) 


d2(a,c) 


-c 

e 


2 

k=0 


c  p  f%  ^  \  <a  1  ^  t *  i 

C(k+l-c)  -(k-c)  ], 


-c  -1 


e  c 


00 

2 

k=l 


l<a<2, 

0<a<l , 


-c 


e 

(a.c) 


2 

k=0 

k/c-1 


[Ik+l-cl^-lk-cT"1], 


c 


1  . 


l<a<2, 

0<a<l , 


and  x^^=|xpsgn(x) . 

Therefore  if  the  resample  size  m  satisfies  m  /n  — ►  c  €  [0,»),  the 

n  n 

averaged-bootstrap  algorithm  works  with  the  appropriate  standardization  given 
in  Theorem  3,  using  appropriate  adjustments  on  the  skewness  and  scale 

parameters  as  in  Theorems  6  and  7,  provided  F  €  (S(a,/3,s)).  When  m^/n  — * 00 , 

the  averaged-bootstrap  distribution  eventually  loses  all  skewness  and  becomes 
symmetric.  In  this  case  the  averaged  bootstrap  works  when  F  €  D^(S(a,0, s)) , 
again  using  appropriate  standardization  and  adjustment  on  the  scale  parameter 
as  in  Theorem  6. 

Figures  1  through  5  illustrate  the  behavior  of  the  distortion  factors  d^ 
and  d2>  as  functions  of  a  €  (0,1)11(1,2)  and  of  c  >  0.  For  fixed  a,  there  is 

essentially  no  distortion  (i.e.,  dj  ~  1  and  dg  ~  1)  when  the  resampling 

proportion  c  is  sufficiently  small  (although  the  rate  of  convergence  may  be 
very  slow,  e.g.,  Figure  3);  this  agrees  intuitively  with  the  case  mn/n  — *  0. 
Distortion  in  the  scale  parameter  can  also  be  nearly  eliminated  by  taking 
mn/n  ~  2/3,  regardless  of  the  value  of  a  €  (1,2)  (see  Figures  1  and  3).  For 
a€(l,2),  the  distortions  in  both  scale  and  skewness  increase  dramatically  as  c 
initially  increases  from  zero;  for  c>l,  these  distortions  remain  fairly  stable 
(see  Figures  3  and  5).  In  general  there  is  less  distortion  as  we  approach  the 


i.o 
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Gaussian  case,  i.e.,  a  f  2  (see  Figures  1,2, 3, 5). 


5.  PROOFS 

5.1.  PROOF  OF  THEOREM  1. 

We  provide  the  proof  only  in  the  case  where  F  has  finite  mean;  when  the 
mean  of  F  does  not  exist  or  equals  ±»,  the  proof  is  essentially  the  same,  in 
fact  simpler  since  cr=  0. 

Athreya  (1987a),  Theorem  A,  shows  that  (1.3)  is  equivalent  to 
C5-1)  {%  (tj)*  j=1 . k)  ^{-PhC^).  j=l, . .  •  ,k) 


for  all  k  £  1,  tj  €  IR,  j=l . k,  where 

.  * 

itT  .  m 

%(t):=E(e  n|*n)  =  [1+ (t)]  n, 

n  n 

with 


m  _ 

:=  r  ^=i  MV  <Yxn» 


(5.2) 


=  (  X  + 
l*Ue  I 


X&)ft(x)dNn(x)  =:  nt£J(t)+Jn>e(t), 


i  tx 

ft(x)  =  e  -1-itx,  for  any  Borel  set  A  C  IR, 

m  1 

(5.3)  N  (A)  :=  &  1I.(a"1(X.-X  )}, 

nv  '  n  j=l  A1  n  v  j  ny J 


and  for  any  x  €  IR,  Nr(x)  :=  Nn([x,«))  if  x  l  0,  and  Nn(x)  :=  Nn((-«,x])  if 
x  <  0.  In  turn  (5.1)  is  equivalent  to 

(5.4)  . k>  S3  (los  >1 . k) 

(see  Theorem  5.5  in  Billingsley  (1968)),  and  by  the  Cramer-Wold  device,  to 
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j 


fjVntV  3s  los  ¥*j> 


for  arbitrary 

We  find  the  necessary  and  sufficient  conditions  for  k=l,  i.e.,  for 
q> 

\/»n(t)  log  <pjj(t)  for  fixed  t  €  1R.  The  conditions  and  the  proof  are 
essentially  the  same  for  any  k. 

In  a  series  of  lemmas  we  find  necessary  and  sufficient  conditions  for  the 
convergence  in  law  of  all  finite  dimensional  distributions  (fdd’s)  of  J( t ) 

XI  t  C 

and  J  (t)  in  (5.2),  and  then  prove  that  all  fdd’s  of  their  sum  converge  in 

n  i  & 

law  if  and  only  if  those  of  both  J(t)  and  J  (t)  converge  in  law. 

We  define  ordered  convergence  in  law,  and  in  probability,  of  doubly 
indexed  random  variables  X  as  follows: 

n  i  & 


itX 


if  -  lim  lim  X  =  X  iff  V  t  €  R,  lim  lim  E  e 


e-O  n-*° 


n.e 


n,e  „  itX 
=  E  e 


6-0  n-*» 


P  -  lim  lim  X  =  X  iff  V  5  >  0,  lim  lim  P( |X  -X|  >  6)  =  0. 


e-O  n-*»  n,fc 


6-0  n-*° 


n,6 


It  can  be  easily  verified  that  many  classical  results  for  r.v. ’s  with  one 
index  are  applicable  here  for  the  ordered  limit,  such  as:  convergence  in 
probability  implies  convergence  in  law,  Slutsky’s  Theorem,  and  the  Central 
Limit  Theorem. 

The  first  two  lemmas  provide  the  necessary  and  sufficient  conditions  for 

the  convergence  in  law  of  the  fdd’s  of  J(t)  as  n  — ►  09  and  then  e  — *  0. 

n )  & 


Lemma  1.  For  each  fixed  t  >£  0, 


P  -  lim  lim 


n.e 


J(t) 


e-O  n-<»  -S 


\x\ie  2 


t^x^dN  (x) 
nv  ' 


=  1. 


Hence  all  fdd’s  of  n  £J(t)  converge  in  law  as  n  — » 00  and  then  e  — »  0  iff 
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<5-5>  U  *?«„(*>  -  r  jj  an2  C*rV2  I(an‘  |XJ-Xnl  *  £> 

o 

converges  In  law,  as  n  — ► 00  and  then  e  — »  0,  and  the  former  limit  equals  -t  /2 
times  the  latter. 

_  -2 

Since  X  converges  a.s.  to  the  finite  mean  of  F,  and  m  a  — r*  0  by  (1.0), 

n  n  n  n-*°  J  v  ' 

s tandard  arguments  enable  us  to  delete  X  in  (5.5).  Therefore  all  fdd’s  of 

n  ' 

-  R  1 

J(t)  converge  in  law  if  and  only  if  X  :=  2  XJ  converges  in  law,  where 

“*  *  ^  |  fe  t  4  H  |  £ 


n,6 


IN  O 

r  &l  X  #  I  fa 
n  n  j  v  1 


-1 

n 


|Xjl 


<  e) 


Note  that  (X^  f°rm  an  infinitesimal  array  of  independent  r.v.’s,  so  we 

can  apply  Central  Limit  Theorem  4.7  in  Araujo  and  Gin6  (1980)  to  obtain  the 
following  characterization. 


Lemma  2.  X  converges  in  law  as  n  — »  »  and  then  e  — *  0  if  and  only  if 

11 1  c 

conditions  (l.b)  and  (l.c)  are  satisfied.  In  this  case  the  characteristic 

2 

function  of  the  limit  Z  is 

<P  2(u)  =  exp{iua2  +  J*  (eiux-l)du(x)}. 

Z 

Corollary.  All  fdd’s  of  J( t)  converge  in  law  if  and  only  if  (l.b)  and  (l.c) 

ll  t  c 

12  2 

are  satisfied,  and  the  limit  equals  -  ^ t  Z  . 

Note  from  (l.b)  that  v  =  0  when  mn/n  is  bounded,  because  the  event 

{  |x|  £  ean,  |x|  >  an  v^nTiiP  }  eventually  becomes  null  for  any  fixed  positive 

2  2 

8,  when  e  is  sufficiently  small.  In  this  case  Z  -a  a.s.,  while  in  general 
2  2 

Z  -a  is  infinitely  divisible  with  Levy  measure  v. 

We  now  find  necessary  and  sufficient  conditions  for  the  convergence  in  law 
of  all  fdd’s  of  J  ft).  It  follows  from  (5.2)  that  the  sample  paths  of 

li  t  w 
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exp {\p  (t)}  (or  exp{J  (t)})  are  characteristic  functions  of  infinitely 
divisible  laws  with  L6vy  measures  the  corresponding  sample  paths  of  Nn(x)  (or 
Nn(x)I( |x|>e)) .  We  show  in  Lemma  3  that  the  convergence  in  law  of  all  fdd’s  of 
Jn  fe(t)  is  equivalent  to  the  convergence  in  law  of  all  fdd’s  of  Nn(x).  Hence 
one  can  focus  attention  on  Nn(x),  which  is  much  simpler  them  ^R(t)  and 
essentially  equivalent  to  a  sum  of  independent  random  variables. 

Lemma  3.  Under  Assumptions  (1.0)  -  (1.2),  the  following  are  equivalent. 

(i)  For  all  k  £  1  and  tj . t^  €  IR, 

S£-  lim  lim  (J  (t.),  j=l . k)  =  (J(t.),  j  =  1 . k}. 

eO  n-*>  n,e  J  J 

(ii)  For  all  k  £  1  and  x^,...,^  €  R\{0), 

{Nn(Xj),  j=l . k)  (N(Xj),  j=l . k}. 

(iii)  Condition  (l.a)  in  Theorem  1  holds. 

Under  any  of  these  conditions,  the  limit  in  (i)  is  given  by 

«C»  if  mn/n  “ 5 *  o. 

J(t)  =  '  c/1  f t(x)dN(x)  if  mn/n  — ►  c  €  (0,»), 

0  otherwise, 

* 

where  N  is  a  Poisson  random  measure  with  intensity  measure  c  *A. 

Notice  from  Lemmas  2  and  3  that  at  least  one  of  J(t)  or  J  (t)  will 

n,e  n,& 

have  a  nonrandom  limit,  for  all  different  choices  of  resample  size  mn» 

Therefore  the  sum  ^  (t)  will  converge  in  law  when  both  of  them  converge. 

Lemma  4:  Under  Assumptions  (1.0)-(1.2),  (5.4)  is  equivalent  to  the  convergence 
in  law  of  all  fdd’s  of  J(t)  and  J  (t),  and  if  (5.4)  holds  then 

n  C  '  '  fl  A  '  ' 
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log  <pu (t)  =  lim  lim  J(t)  +  lim  lim  J  (t). 
H  e-O  nn»  n'e  e-O  n-*»  n,e 


Proof  of  Lemma  1.  Fix  6  >  0.  Then,  for  each  n  and  e  <  35/|t|  we  have 

,J(t) 


n,£ 


-J I  I,  (x) 

|x|£e  2  nv  ’ 


-  1 


>  6 


=  p[  U{xHe  (eltx-l-itx  ♦  it2x2)dKn(x)|  >  5  X|xUe  |t2x2dMn(x)  ] 

*  PC  J'|x|$t  g;l tx |3dNn(x)  >  6  SMie  |t2x2dNn(x)  ] 

*  pC  |t2(||t|e-6).  X  x2dNri(x)  >  0  ]  =  0, 


|x|£e 


because  |elx  -  1  -  ix  +  x2/2  |  $  |x|2/6  for  all  x  €  R. 

Proof  of  Lemma  2.  We  first  verify  that  (X^  ,  l£j£n)  form  an  infinitesimal 

ri)  £ 

array.  Indeed  for  fixed  e  and  6  >  0,  we  have  by  (1.0), 


P(X^  t  >  6)  =  Pfa”1  |x|  .<  e.  ^  a^2  >  8)  i  P(|x|  >  a ^  0. 


n,6  ~  n  -n  -  n  n'  n-*» 

Now  by  Theorem  4.7  in  Araujo  and  Gin6  (1980),  2I?_1X^  converges  in  law 

j— i  n  t  & 

(as  n  — ►  80  and  then  e  — *  0)  if  and  only  if  the  following  three  conditions  are 
satisfied: 

(i)  there  exists  a  LAvy  measure  i>  such  that  for  all  0  £  6  <  u  with 
u{6)  =  0  =  u{u), 


lim  lim  n  P (6  <  X1  i  u)  =  v(5, u]; 

6-0  n-*9  n,e 

o 

(ii)  for  all  6  >  0  such  that  v{6}  =  0,  there  exists  cr^  ^  0  such  that 


lim  lim  n  E[X* 
e-O  n-*»  n,£ 


I(X 


n,e 


< 


5)1 


2 

a6i 


(iii)  there  exists  a  €  IR  such  that 
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lim  lim  11m  n  VarCX1  I(X*  <  6)]  =  a2; 

6-0  e-0  n-*»  n,e  n*e 

in  this  case  the  characteristic  function  of  the  limit  is 

exp{iucr2-  |u2a2  +  J,Q(eiux-'.)du(x)} , 

2  2  2 

where  a  =  lim^jQ  a g,  which  exists  since  a g  is  monotonic  (decreasing)  6  i  0. 

Conditions  (i)  and  (ii)  are  seen  to  be  the  same  as  (l.b)  and  (l.c)  in 
Theorem  1,  and  (iii)  is  implied  by  (ii)  since 

n  <  5»2  i  <  5>]  S3  rr*  0  . 

6.-0  5i0 

2 

Hence  a  =0,  i.e.,  the  Gaussian  component  in  the  limit  is  degenerate,  and  the 

limit  has  ch.f.  as  shown.  □ 

Proof  of  Lemma  3.  (ii!  is  equivalent  to  fiiil;  We  will  show  that  (iii)  is 
equivalent  to 

(5.6)  (Nj^Cyj,00),  NnC-00*-}^])  converges  in  law  for  all  ylty2  >  0, 

since  the  multivariate  argument  can  be  similarly  carried  out  to  achieve  (ii). 

Again  by  a  standard  argument,  we  can  drop  X  in  N  ,  so  for  any  y.  ,y0  >  0, 

n  n  l  ^ 

Uj.Ug  €  R, 

(5.7)  Eexp{i(u1Nn[y1,»)+u2Nn(-«,-y2])} 

1,1  1 

=  Eexp{i  JL  ^=1[u1I(a;1Xj  l  yx)  +  u2  Ka^X.  <  -y2)]} 
m  1 

=  C  Eexp{i  ^  [UjIfa^X  *  +  u2  I(a^X  i  -y2)]}  ]n 

=  Cl  +  J  (  CexPdUjin/n)  -ljPfa"^)  +  [exp(iu2mn/n)  -l]P(a^X^-y2)}  ]n. 
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Hence  (5.6)  is  equivalent  to 

n[exp(iumn/n)  -l]P(a~*X  £  y)  and  n[exp(iumn/n)  -l]P(an*X  <  -y) 


converging  for  all  y  >  0,  u  €  IR,  which  is  equivalent  to  (l.a).  Therefore  (5.6) 
is  equivalent  to  ( i i i ) .  Moreover  (5.7)  converges  to 

*  exp{i[u1X[y1.®)  +  u2\(-«,-y2]]}  if  nwTi  — ►  0, 

■  exp{c  ^expfiUjCj-lp^yj,®)  +  c  1[exp(iu2c)-l]X(-«,-y2]}  if  mn/n  — >  c€(0,«>)  , 

1  otherwise. 


(il  implies  fii):  For  simplicity  of  notation  we  prove  (ii)  is  implied  by 
the  convergence  in  law  of  all  fdd’s  of  ^  (t). 

From  (5.2)  we  can  write,  with  g^(u)  =  e*tU-l, 

V*)  =  C.  ft(x)dNn(x) 


(5.8) 
where 

(5.9) 


=  S0  St(u)du]  ^(x)  +  J^C-itJ^  gt(u)du]  dNn(x) 

=  gt(x)Nn[x,»)dx  -  itJ'^<ogt(x)Nn(-®,x]dx 
=  -t 2jT4a  eitxMn(x)dx  =:  -t2^(t), 


M  (x)  = 
nv  ' 


'  C  vu’“)du’ 

.  ■£»  Nn(-»,u]du, 


x  >  0, 

x  <  0. 


Thus  yp ■  ( t)  is  the  Fourier  transform  of  Mn(x),  which  is  piecewise  linear, 
continuous  with  compact  support. 

It  is  clear  from  (5.2)  that  expf^ft)}  is  a  (random)  characteristic 
function  of  an  infinitely  divisible  law  with  (random)  L6vy  measure  Nn>  and  (i) 
implies  the  fdd’s  of  yp  ( t)  converge  in  law  to  the  fdd’s  of  log  ^jj(t)  =’■  '/'(t). 
If  the  ^n’s  were  nonrandom  or  more  generally  if  the  convergence  were  a.s. 
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instead  of  in  law,  the  convergence  of  their  L6vy  measures  in  (ii)  would  follow 
immediately  from  classical  results. 

In  order  to  establish  (ii)  in  the  current  case,  we  apply  the  embedding 
theorem  of  Skorohod  (1956)  to  the  ^’s  and  'p,  whose  sample  functions  are  in 
C(IR).  In  this  special  case  of  processes  with  continuous  paths,  this  implies 

(V  /v 

that  there  exist  stochastic  processes  and  'Kt)  on  some  probability 

space  (n,?,P),  such  that:  the  fdd’s  of  'p  and  'p  are  the  same,  and  so  are  those 
of  and  for  each  n,  and  the  sample  paths  of  ^(t)  and  'P(t)  are  with 
probability  one  in  C(IR)  (therefore  their  corresponding  induced  measures  on  C(IR) 
are  the  same) ,  and 

(5.10)  *n(t)  for  all  t. 

Therefore  there  exists  0^  C  fi  with  P(12j)  =  0,  such  that  for  all  w  € 

(5.10)  holds. 

We  will  show  that  each  >1/  has  the  same  form  as  'p  ,  so  that  exp {*p  }  is  the 

rn  n  nJ 

characteristic  function  of  an  infinitely  divisible  law  with  Levy  measure 

having  the  same  fdd’s  as  Nn.  Then  (5.10)  will  imply  a.s.  convergence  of  Nn> 

therefore  the  convergence  in  law  of  all  fdd’s  of  N  . 

Since  ^  and  have  the  same  induced  measures  on  C(R),  in  order  to 

establish  that  \p  has  the  same  form  as  'Pn>  determined  by  (5.8)  and  (5.9),  it 

suffices  to  show  that  the  inverse  of  the  map  that  defines  \p  or  from  N  via 

r  n  n  n 

Mr  is  measurable.  We  now  proceed  to  discuss  this. 

We  first  formalize  the  two  maps  from  N  to  M  ,  and  then  from  M  to  \p  . 

r  n  n  n  n 

Let  M  (or  Jf)  be  the  set  of  all  discrete  signed  measures  (or  discrete  measures) 
on  IR  with  finitely  many  atoms  with  the  topology  of  weak  convergence.  Let  ?  be 
the  set  of  all  functions  that  are  piecewise  linear  with  finitely  many  segments, 
continuous,  with  compact  support.  2?  is  a  measurable  subset  of  C(IR)  with  the 
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uniform  topology.  Define  T^:  A  — *  9  by 

£  u[y.w)dy.  x  i  o, 

(V)(x)  =  J 

L  u(-».y]dy.  x  <  0. 

for  all  v  €  A.  is  one-to-one  and  onto  and  its  inverse  9  — >  i  is 

defined  by  (T^fXa.b]  =  (Bf )(b)-(Bf )(a) .  where  B.B^:  9  — »>  K(K)  are  given  by 

(Bf)(x)  =  lim  (B  f)(x).  (B  f)(x)  =  £  [f (x+r,)-f(x)]. 

77X0  v  71  v 

B  is  a  measurable  map  since  B  is  continuous  for  each  fixed  tj  >  0.  Therefore 

17 

is  well  defined  and  measurable. 

The  relationship  (5.8)  is  the  Fourier  transform  Tg:  Cq(IR)  — *  C(IR) ,  where 
C(K)  (or  Cq(IR))  is  the  set  of  all  continuous  functions  on  IR  (with  compact 
support).  T2(C0(R)}  =  T?(ir= jCgt-m .  m]  }  =  ^.^{CQC-m.m]}  is  a  Borel  subset  of 
C(R),  since  the  sets  C^f-ra.m]  of  all  continuous  functions  with  support  on 
[-m,m]  are  closed  subsets  of  C(IR),  and  T2  is  continuous  on  each  CQ[-m,m]  as 


,,T2(fl>"T2(f2)H«  i  «flf1(x)-f2(x)|dx  S  2m  llfj-f^. 


Hence  we  can  define  the  inverse  map  T2*:  T2(Cq(R)}  — »  Cq(IR)  by 

1  00  ^  10  0  ^ 

(T2Af)(x)  =  lim  e  ltx  f(t)exp{-  |oV}dt  =:  lim  (A  f)(x). 

CT-O  CT-O 

for  each  f  €  T2(Cq(R)}  (see  e.g. ,  Rudin  (1966),  §  9.7).  T2*  is  Borel 

measurable  since  each  map  is  continuous  on  T2(Cq(IR)}: 


"V?i>  -  A„Cf 2> i  lfr'2"»  h  C  «*<-  i°2t2>dt  -  5 w  "Vy»  • 


Therefore  we  have 
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M  C  J  fPTt  »  C  C0(K)  fTTf  T2{C0(E)>  C  C(E). 

X1  *2 

Then  B  =  T^T^A)}  =  (T”1  )—1  { (T"1 )— 1  (^) >  is  a  measurable  subset  of  C(IR).  Since 

the  paths  of  all  ^*’s  are  in  B  with  probability  one  by  (5.8)  and  (5.9),  we  have 
n 


i = p  .  *>> = ?  •  cv>. 


which  implies  '/^(•.w)  belongs  to  B  a.s.  (P).  Thus  there  exists  Aq  C  17  with 
P(Aq)  =  0,  such  that  for  all  u>  €  A\Aq,  we  have 

Nn(-,u)  :=  ^(T-1^.,.)]}  €  M ,  n  =  1.2,... 

i.e. , 

*n(t)  =  Xft(x)^n(*)»  n=l ,2, . . . 


Therefore  for  all  w  €  A\(AjUAq)  and  all  t,  we  have 
(5.11)  exp{XftdNn)  =  exp$n(t)}  ^  exp{^(t)}. 

»v  /x/ 

Applying  the  classical  result,  the  Levy  measures  of  exp^^t)}  must 

converge  a.s.  (P)  to  the  Levy  measure  N  of  exp{^(t)}. 

Now  notice  that  N  and  N  have  the  same  fdd’s  because  for  all  Borel 
n  n 

subsets  A  of  K(IR) 

P(Nn  e  A)  =  P0/£  efT"1)"1^"1)"1^)})  =  P(**  efT"1)"1^”1)"1^)})  =  P(Ne  A). 
Therefore  for  all  x^ . x^  €  IR\{0), 


<v 

(N^),...,^^))  I  (Nn(Xl) . \(\))  (Nfxj),...,^)) 


A/  ^ 

which,  implies  (iii)  by  the  previous  argument.  Hence  N  =  N  for  some  Poisson 

A/ 

random  measure  N  with  intensity  measure  A.  Since  N  is  a  random  Levy  measure, 
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x^dN(x)  <  »  a.s.  (P)  and  thus  J^x^dN(x)  <  »  a.s.  (P).  Therefore  the  form 

1  2 

of  its  ch.f.  implies  that  J*_jX  dX(x)  <  w,  and  hence  X  is  a  Levy  measure. 

(ii!  implies  fils  We  write 


+  J  „(t). 
n,e,Mv  ' 


For  the  first  term  we  have 


n  m 


=  2 1 1 1  2  —a  1|X.-X  |  I(a-i|X,-X  |  l  M) 

1  1  n  n'j  n'  vn'j  n,fc  ' 


_“1  Iv  V 


n  m  1  .  n 

=  2 1 1 1  2  —  a^jxj  I(a“A|X.|  \  M)  =:  2|t|  2  f£ 

1  1  j  j  n  n  1  j 1  v  n  1  j 1  u  1  1  1  j  n. 


M 


\n 
n, 

infinitesimal  array  as  n  -»  »  and  then  M  -»  09 ,  since  for  any  5  >  0, 


because  by  a  standard  argument  we  can  drop  the  sample  mean.  {W^  *s  311 


P(*ifM>  5)  =  P(|X|  *  anM.  |X|  *  6ann/mn)  £  P(|X|  *  anM)  ^  0. 


Therefore,  applying  the  degenerate  version  of  the  C.L.T.  (Loeve  (1977), 
p.  329),  we  have 


(*) 


n 

P  -  lim  lim  2 

M-s»  n-*»  j=l 


=  0 


if  and  only  if  the  following  three  conditions  are  satisfied: 


(a)  lim  lim  n  ECW1  u  I(WX  „  i  6)]  =  0  for  some  6  >  0, 

ii  ^  ^  n,M  '  n,w 

M N»  nH» 

(b)  lim  lim  n  PfW1  M  >  6)  =  0  for  all  6  >  0, 

n n*n 

(c)  lim  lim  n  Varrw1  u  IfW1  u  £  6)]  =  0  for  some  6  >  0. 

M-*>  n-*»  n*M  n*M 
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It  can  easily  be  verified  that  (a)  and  (b)  are  equivalent  to  assumptions 
(1.1)-(1.2),  and  (c)  is  implied  by  (a)  due  to: 


VarK.H  ^  <  EC<,K  Wn.M  *  ‘  S  *  5»' 


It  follows  that  (*)  holds.  Hence  the  fdd’s  of  J  (t)  converge  in  law  (as  n-**> 

and  then  e  -»  0)  if  and  only  if  those  of  J  M(t)  do  (as  n-*»,  then  e  -*  0,  and 

n  1 6 1  ri 

then  M  -»«),  and  their  limits  are  the  same. 

Now  ft(x)  is  continuous  and  bounded  on  [-M,M]  for  each  fixed  M  >  0,  hence 

the  convergence  in  law  of  all  fdd’s  of  N  implies  that  of  J  ,  u  as  n  -*  ®, 

n  n.e.w 

e  -*  0,  and  then  M  -»  ®: 

J  M  -4?  J*  f  .dN  J*  f  dN  J(t).  0 

n,e,M  n-J®  t^|x|<M  z  ^  |X|<M  t  M 


Proof  of  Lemma  4.  (5.4)  implies  the  convergence  in  law  of  all  fdd’s  of  N,  and 
hence  of  £  (Lemma  3) .  Notice 


V*) 1  V‘>  ■  . . 

MS 


+  X  )f  tdN  =: 

e  |x|>e 


n.e 


J(t)  + 


w*)- 


where  the  fdd’s  of  ^  (t)  converge  a.s.  (P)  to  those  of  'P(t)  (see  Lemma  3). 

/V  A/ 

Since  N  converges  a.s.  (P)  to  N,  the  Levy  measure  of  exp{'//(t)}f  J  (t)  will 
n  iii  6 

converge  a.s.  (P)  as  n  — *  ®  and  then  e  — ►  0  because  Nr(x)I( |x|  >  e)  is  again  a 
Levy  measure  for  fixed  e  >  0,  and 


lim  lim  J  (t)  =  lim  X  f .dN  =  ff  dN  a.s.  (P). 
e-*0  n-s®  n,e  e-O  |x|>e 


Therefore  n  feJ(t)  =  $n(t)  -  £(t)  converges  a.s.  (P)  as  n  — »  ®  and  then 

e  — ¥  0  to  'p(t)  -  J'f^.dN.  (Note  'P(t)  may  have  a  random  Gaussian  component 
-  Jjt^a^(w)  which  distinguishes  'P(t)  from  Xf^dN).  Furthermore, 


J(t)  =  X  f  dN  I  X  f  .dN  =  J(t) 
n,e  v  1  it.  tn  i  i .  t  n  n.e  v  * 
|x|Se  |x|Se 
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since  and  Nr  have  the  same  fdd’s  for  all  x  €  IR\{0},  and  ft(0)  =  0. 

Therefore  the  fdd’s  of  J(t)  must  converge  in  law. 

n  •  £ 

Conversely,  if  the  fdd’s  of  J(t)  and  J  (t)  converge  in  law,  then  so 

do  those  of  their  sum  J(t)  +  J  (t),  since:  In  each  of  the  two  mutually 

n.fe  n,£v 

exclusive  cases  (i)  m  /n  — »  c  €  [0,°°),  and  (ii)  otherwise,  either  J(t)  or 
v  y  n  n,e 

J  (t)  converges  in  law  to  a  nonrandom  limit,  and  therefore  converges  in 
n  i  £ 

probability,  thus  ensuring  the  convergence  in  law  of  the  sum  by  Slutsky’s 
Theorem.  □ 


5.2.  PROOFS  OF  THEOREMS  2-8. 

y> 

Proof  of  Theorem  2.  Since  for  each  fixed  x  €  IR,  (x)  —  —  >  H(x) ,  and 

n 

0  £  ,  H  £  1,  it  follows  that  EHn(x)  =  EjH^  (x)  —  >  G(x) .  On  the  other 

n  n 

hand,  by  the  inequality  of  Hoeffding  (1963),  for  all  e  >  0, 

1  € 

P(|Hn(x)  -  EHn(x)|  >  e)  =  P(  j-  ^  (x)  -  E^  (x)]  >  e) 

£  2exp{-2£^eP'}  £  1/n^, 

for  sufficiently  large  n.  The  Borel-Cantelli  Lemma  implies  a.s.  convergence  of 
Hn>  and  uniform  convergence  follows  in  the  standard  way.  0 


Proof  of  Theorem  3.  Put  x(t)  =  l--F(t)+F(-t)  for  t  >  0. 
(i)  a  =  2.  We  have  (Feller  (1966),  p.  545) 

(5* 12)  Urn  - =^=0' 

t-*>  Epri(|X|  *  t)]  a 

Case  (a).:  mn/n  — *  c  €  [0,»).  (l.a)-(l.c)  imply 

mn^(any)  XCy’°°)  + 
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VvV2  I(|x|  $  ean)]  =  ^ o 2  ■ 

e-O 

These  imply  \  =  0  by  (5.12),  and  m  a  ^(a^  — *  ex2  (>0). 

Case  f  b) :  m  /n  — *  00 .  Put  b  =  a  vWm  .  Then  b„  — ► 09  by  (1.0).  Also 
- *— L  n  nnn  n  v  ' 

(l.a)-(l.c)  imply  for  all  6  >  0,  n\(bnv£f  )  — *  v(5,*>),  and 

a?=lim  m  a-2EX2I(|x|  ^  b  V<5  )=lim  m  a“2L(b  v£)  =lim  nb'^fb^)  =lim  nb“2L(b  ). 
5  n-*>  n  n  n  1  n-»»  n  n  n  n-*>  n  n  n-*»  n  n 

Again  using  (5.12),  we  have  v  =  0,  and  rib  ^Lfb  )  a2  £  0.  In  either  case, 

H  is  Gaussian. 

(ii)  0  <  a  <  2.  Case  fa):  m^/n  — »  c  €  [0,®).  (l.a)  implies 

m  (ya  )  ^(a  y)  converges  for  all  y  >  0;  thus  m  a  °TL(a  )  — ►  constant  Q  0). 
n  n'  n  nnn 

Case  fb) :  m  /n  — ►  ».  (l.b)  implies  n  x(b  v^f)  ~  (c1+cQ)n(b  v?f)  ^(b  V5) 

n  n  x  £»  n  n 

converges  for  all  5  >  0.  □ 

Proof  of  Theorem  4.  By  Athreya  (1987a), 

<PG(t)  =  Eexp(J^w  ft(x)dN(x)}, 

where  N(A)  =  N'(At),  N'  is  a  Poisson  r.m.  with  intensity  A^,  and  r  is  the  last 
jump  of  N 1 .  Hence 

<<>G(t)  =  exp{f t(l)}  Eexp^^  1^ft(x)dN(x)}  =:  exp(ft(l)}  £(t). 

It  can  be  shown  that,  given  t,  N  is  a  Poisson  r.m.  with  intensity  measure 
Xa(xT)  =  t  “x  (x) ,  and  that  t  has  max-stable  distribution  for  which 

Eexp(uT-a)  =  c^a  exp{(u-c1)y  “}  y  a  *dy  =  [l-c^u]  1 

if  Re(u)  <  Cj.  Hence 

*(t)  =  E  [  E  (exp{  X  f  (x)dN(x)}  |  r  }] 

(-«,1) 
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=  E  exp{T~aX^  [exp{f t(x)}-l]  dXa(x)} 

=  [  1  -  c^X^  [exp{f t(x)}-l]  dXa(x)  3"1, 
since  the  real  part  of  the  last  integrand  is  £  0  <  c^.  0 

Proof  of  Theorem  5.  In  this  case,  Athreya  (1986a)  showed 

<PH(t)  =  exp{/^ft(x)dN(x)}. 

where  N  is  a  Poisson  r.m.  with  intensity  Therefore 


</>G(t)  =  exp{  SZa  Cexp{ft(x)}-l]dAa(x)  }. 

When  1  <  a  <  2,  we  find 

-log*  (t)  =  a[ClX“  +  c/j(l  -  exp{f  (x)})  — 

M 

=  a|t|a  [  CjX*  {1  -  exp[cosu  -1  +i(sinu-  u)sgn(t)]  }-~y^ 

u 

+  c2X0  {1  -  exp[cosu  -l-i(sinu  -  u)sgn(t)]}-r^  ] 

u ' 

■  M“a(«i«2)CJ  '  i  «n(t)3. 

where 

a  =  a  Xq  [1  -  exp(cosu-  l)cos(sinu  -  u)]  ■  •  , 

u 

b  =  a  Xq  exp(cosu  -  l)sin(sinu  -  u)  -y— ■  . 

u 

Integrating  by  parts  and  then  using  the  series  expansion  of  exp{e*u),  we  obtain 

oo  (x 

T(2-a)  ,  n  car  t.ocn^r  -o  k 

a  -  ib  ”  (a-l)e  ^  ~2  cos  2  “  exp^2  (k+1)! 


-  1]  }• 
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hence  b/a  =  tg(cnr/2)d2(a) .  Since  a(cj+c2)  =  sdj(a),  /3  =  (ci-c2)/^ci+c2)  • 
s  =  (Cj+c2)r(2-a)  |cos(air/2)  |/(a-l) ,  it  follows  that 

<PG(t)  =  exp{  -sdj(a)  1 1  |a[l  -  i/3d2(a)tg(arr/2)sgn(t)]  }. 

The  proof  for  0  <  a  <  1  is  similar  with  a  -  ib  =  e  ^^^^(l-ajdjfa) .  0 


Proof  of  Theorem  6.  If  a  =  Tn  d„ fa  11  ,  then  we  have 

-  n  i  n 

A 

a"1  n1/a  =  exp{  (aan)-1(an-a)log  n  }  dj  n(an)  ^  d ^/a{a), 

—  1  /a 

for  dj  (a)  is  continuous  over  (0,1).  Therefore,  for  A  =  [x,«),  we  have 

P[N  (A)  =  k]  =  P[a-1 (X  lAl-X  )  l  x.  a-1(X  .  -  X  )  <  x] 

L  nv  J  J  L  n  v  n,n-k+l  n/  n  v  n,n-k  nJ  J 

nr  1/ CL  “1  /CL***  T7  \  v  1  /CL  “1/ tt/y  y  \  s  *1 

-  pCan  »  »  (Xn.n-k+l  -  Xn>  >  *•  an  ”  "  (Xn,n-k  “  Xn>  <  x] 

Sp("(*)  -  k). 

by  Slutsky’s  Theorem,  where  N  is  a  Poisson  r.m.  with  intensity  X^d^Ca). 
Therefore  it  can  be  shown  (see  Proposition  1,  Athreya  (1987a))  that 


P  -  lim  lim  J(t)  =0,  P  -  lim  lim  J  M(t)  =  0. 


e-O  nn» 


M-*»  n-J» 


*n(t)  ^|+/ft(x)dN(x), 


and  similarly  for  their  fdd’s. 

Applying  the  proof  of  Theorem  5,  we  find 

<PG(  t)  =  exp{  -sd^  (a)d  j  (a)  1 1  |a[  1  -  i/3tg(a7r/2)sgn(t)]  }.  0 

Proof  of  Theorem  7.  Since  X  =  X+-  X  is  in  D(S(a,/3, s)) ,  it  follows  that 
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Z  =  b  X  -  b  X  is  in  D(S(a,/3z>  s^)) ,  where 


Note  that 


pz  =  (Clb+-c2b  )/(Clb++c2b")  =  p/ d2. 
sz  =  sfc^+Cgb  )/(c1+c2)  =  s(l-/32). 

%  <*>  -  <  1  +  k  3=1  ftta;l(ZrZn)]  }» 


converges  in  law  if  and  only  if 


(5.14) 


,  f Ta-1(Z.-Z  )] 
j=l  tL  n  v  j  n/J 


converges  in  law.  Notice  {Zj}j_^  are  no  longer  iid  due  to  the  estimates  and 

A 

Pn>  which  depend  on  the  entire  sample  9^.  Hence  we  must  show  that  (5.14) 
asymptotically  behaves  as  when  a  and  P  are  known;  then  we  can  handle  the 
problem  similarly  as  in  Theorem  5.  The  arguments  are  as  follows. 

Let  Zqj  =  b+Xj  -  b  Xj,  Z^  =  b+X*  -  b  X^.  Then  using 
I f‘t:(x)— £ t(y)  |  £  2|t||x-y|,  and  writing  b*  for  b*(an,0n) ,  we  have 

l  -  3=i  l 

i  2|t| 

i  Zl'la^d^l  ^IX^I  +  |b--b- 1  |x--ri> 

(5.  IS)  $  [a;1  a17"]  „-1/a{  |S;-b+ 1  ^=1x;  ♦  l^-b'  I^X'}  4 1 1 1  . 

The  first  factor  converges  in  probability  to  [dj(a)(l-/32)]  using  a  similar 
argument  as  in  Theorem  6. 

When  0<a<l,  the  terms  |bn~b  |,  |bn~b  |  converge  in  probability  to  0 

A  A  | 

because  a^,  p are  consistent  and  b  ( a,p )  are  continuous;  and  the  terms 


34 


jXj  converge  in  law  since  X+,X  also  belong  to  a  domain  of 
attraction  (with  the  same  index  a). 

When  l<a<2,  the  terms  n  converge  a.s.  by  the  SLLN;  and 

A*|.  1— 1/ct  A—  —  * 

the  terms  n  [b^-b  l»  n  Ib^-b  I  converge  to  zero  in  probability  since 

lyV - 1— {  [1-/3-P  -d,(S  )](P  -p)  -  P(l-P)[d2(S  )-d2(a)]/d2(a)  }, 

'  d2<an> 

b'-b"  =  - i—  {  -[1+P+P  -d2(S  )](P  -P)  +  P(l+P)[d2(Sn)-d2(a)]/d2(a)  }, 

d2<“n> 

A  A  1  —  1/CC 

an>Pn  are  consistent  estimates  of  a,p  with  rate  n  ,  and  is 

differentiable  on  (1,2). 

Therefore  (5.15)  converges  to  zero  in  probability,  and  (5.14)  converges  in 
law  if  and  only  if  f t[a”1(Z0J-Z()n)]  converges,  and  they  have  the  same 

limits.  We  can  now  use  arguments  similar  to  those  in  Theorems  5  and  6  on 
(Zqj}j_i  to  establish  the  result.  0 

Proof  of  Theorem  8.  Suffiencv  of  both  fil  and  fill:  Assume  F  €  D(S(a,/3,s)) , 

or  D  (S(a,0,s))  when  m  /n  — *  0. 
p  "  n 

Case  (11:  mn/n  — *  c  €  [0,«°).  Then  by  Theorem  3,  (1.3)  implies  that 

lim  m  a  T.(a  )  =  d  for  some  d  £  0,  and  G  =  H,  which  is  normal  with  a  =d(2,c) 
n-*° 

when  a  =  2;  while  when  1  ^  a  €  (0,2), 


a 2  =  lim  lim  m  a “2  E[X2I( |X|  £  a  e)] 
e-O  n-*°  n  n  n 


=  lim  lim  m  a  (a  e)  ^(a  e)  =  lim  e  lim  m  a  ^(a  )  =  0. 
-jmnnvn-'  v  n  '  _  ^  n  n  v  n' 

e-O  n-*°  e-O  n-*° 

Also  by  (l.a), 

X[y,«)  =  lim  mnP(X  >  yaR)  =  lim  c1mn(yan)"aL(yan)  =  Cjdy  a. 


n-*® 


n-#» 


and  similarly  \(-«,-y]  =  C2dy  a,  so  that  X  =  dXa.  Hence  by  (1.5), 
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*G(t)  =  1 


exP  {<*/_«  ft(x)d^a(x)> 


-lr“ 


if  m  /n 
n 


0. 


I  exp{dc  S_„  [exp{cft(x)}-l]dXa(x)} 


if  mn/n  — »  c  €  (0,«); 


=  « 


if  m  /n  — *  0, 
n 


exp{-ds|t[  [1  -  i/3tg(cnr/2)sgn(t)]} 
exp{-dsdj(a,c)  |t  |a[l  -  i/3d2(a,c)tg(a7r/2)sgn(t)]}  if  n>n/n — *c  €  (0,»); 


using  arguments  similar  to  those  in  Theorem  5.  In  fact  Theorem  5  is  a  special 

case  of  Theorem  8  with  c  =  1,  d  =  1,  d^(a)  =  dj(a,l),  and  d2(a)  =  dgfa.l). 

Case  (2) i  m  /n  — ►  ».  Then  bv  Theorem  3.  lim  nb  ^fb  )  =  d  for  some  d  £  0, 

- —  n  _ _  n  n7 

n-»“> 

where  b  =  a  v4i/m  .  Moreover  G  =  H  is  normal  with  a  =  d(2,c)  if  a  =  2,  while 
n  n  n  v  1 

when  1  ^  a  €  (0,2)  we  have  from  (l.b)  (with  0<6<u  such  that  i>{5)  =  0  =  u{u}). 


u(6,u]  =  lim  lim  n  P(b  v^5  <  |x|  <,  b  v£T  ) 
eO  n-*>  n 


=  lim  (c.,+c0)n[(b  vS)  ^(b  vG)  -  (b  \G)  ^(b  vG)]  =  (c1+c0)d(6  a^2-u  CC^2). 
^  l  /  n  n  n  n  i  ^ 

n-joo 

Condition  (l.c)  can  be  written  as 


o2  =  lim  lim  m  a  ^{X^dxl  $  bvfr)} 
t-O  n-*»  n  n  n 


=  lim  nb'2  (bnvG)2-aL(bnvG)  =  51_a/2  lim  nbn“L(bn)  =  d61_a/2  ^  0. 


n-*» 


n-*° 


Therefore  by  (1.5),  <pg(t)  =  Eexp{-  'At^N} ,  where 

«Pw(t)  =  exp{  d(Cl+c2)|  (eitX-l)  }* 

which  is  the  ch.f.  of  S(a/2,l,ds).  Hence  it  follows  that 


<Pq( t)  =  exp  (-|t|asd/[2a/,2cos(cnr/4)]}. 


Necessity  of  fill:  mn/n  — *  0.  If  a  =  2,  it  follows  from  the  expression  of  <Pq 
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in  (1.5)  that  X  =  0,  i.e.,  for  y  >  0, 


m  P(a  lX  >  y)  -*  0.  m  P(a  *X  <  -y)  — i ►  0. 
n  v  n  JJ  n  v  n  JJ 


Also  there  exists  a  such  that 


o'  =  lim  lim  m  a  2  E{X2  I(|X|  *  ea  )}. 
e-O  n  n  n 


Hence  F  €  Dp(S(2,0,s))  (Ibragimov  and  Linnik  (1971),  Eq.  2.6.12-2.6.13,  and 
Feller  (1971),  p.  555). 

2 

If  1  a  €  (0,2),  then  a  =0,  and  X  is  the  Levy  measure  of  an  a-stable 
law,  which  means  there  exist  constants  c^,  ^  1  0,  such  that 


X[y,»)  =  Cjy"®,  X(-«,-y]  =  c^y~ 


Thus  by  condition  (l.a), 


lim  mnP(a~XX  >  y)  =  lim  ■  P(a“1X  <  -y)  =  c2y“ 

n-*»  nn« 


and  by  (l.c) , 


lim  lim  m  a  2  E[X2  I(|x|  $  ea  )]  =  a2  =  0. 
e-O  n-#»  n  n  n 


Again  these  imply  F  €  Dp(S(a, *,*))  (see  Ibragimov  and  Linnik  (1971),  Eqs. 
(2.6.3)-(2.6.5)  and  Feller  (1971)). 

Necessity  of  f i 1 :  Case  fll:  mn/n  — *  c  €  (0,»).  The  argument  for  a  =  2  is  as 
in  the  necessity  of  (ii).  Now  if  G  is  a-stable  with  1  jt  a  €  (0,2),  then  by 


(1.5), 


'/'(t)  :=  log  fG(t)  =  -  |t2a2  +  c"1/^  [exp{cf t(x)}  -  l]dX(x) 
=  -  |t|as{l  -  i/3tg(a7r/2)sgn(t)}. 


First  we  show  that  a  -  0.  For  t  >  0,  we  have 


(5.16)  h  t2  aa2  +  c  *t  (1  -  exp[c(costx-l)]  cos[c(sintx-tx)]}  dX(x)  =  s. 
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Since  1-  exp[c(costx-l)]cos[c(sintx-tx)]  £  0  for  all  t,x  €  1R,  and  s  €  (0,»), 

o  o-a 

letting  t  — in  (5.16)  implies  a  =  0,  for  t  — ►  «.  Now  exp^t)}  is  the 

1/a 

ch.f.  of  an  a-stable  law,  therefore  it  satisfies  n^(t)  =  \p(n  t),  i.e., 
n  J*Z»texP{cft(x^  "  ^  =  «£«[exP{cft(x^  -  !]  dX(n~1/ox), 

1  /rt 

where  X(x)  =  X[x,«>)  if  x  >  0,  X(-«\x]  if  x  <  0.  Therefore  n  X(n  /  x)  =  X(x). 
By  the  monotonicity  of  X,  we  can  conclude  that  X(y)  is  proportional  to  y  , 
-«0<p<+<»,  (Feller  (1971),  VIII  8,  Lemma  3),  and  p  has  to  be  -a.  Therefore 
F  €  D(S(a, *,*))  by  using  similar  arguments  as  in  the  necessity  of  (ii). 

Case  f2~):  mn/n  — *  “•  Since  G  is  stable,  and  (1.5)  implies  that  it  must 
be  symmetric,  we  have 

E  exp{-  4  t2(cr2+W)}  =  exp{  -s  |t|a), 

M  O 

2 

which  implies  a  +  W  ~  S(a/2,l,Sj),  since  the  left  hand  side  is  the  Laplace 
2  2 

transform  of  a  +  W  at  t  /2.  Therefore 

Eexp(it(a2  +  W)}  =  exp{itcr2  +  Jq  (e*tx-l)du(x)}  =  exp{-Sj | t |a) . 

2 

Hence  a  =  0  and  i>  is  the  Levy  measure  corresponding  to  S(a/2,l,Sj),  i.e., 

lim  n  P(|X|  >b^)=  c (v5)““, 
n-*»  n 

and  with  a  similar  argument  as  in  Case  (1), 

0  =  a2  =  lim  lim  nb"2  E[X2  I( |x|  ^  b  v£)]. 

6-0  n-*»  n  n 


Hence  F|^|  €  D(S(a, *,♦)),  which  implies  F  €  D(S(a,»,*))  since 


F  €  D(S(a, » ,  •))  iff 


t2x(t) 


2-a 


E[X2I(|xUt)]  tH»  a 


(see  Ibragimov  and  Linnik  (1971),  Theorem  2.6.3). 
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